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ABSTRACT 

We report the discovery of 70.5 ms pulsations from the X-ray source AX J1838. 0—0655 using the 
Rossi X-ray Timing Explorer (RXTE). PSR J1838— 0655 is a rotation-powered pulsar with spin-down 
luminosity E = 5.5 x 10 36 ergs s _1 , characteristic age r c = P/2P = 23 kyr, and surface dipole magnetic 
field strength B s = 1.9 x 10 12 G. It coincides with an unresolved INTEGRAL source and the extended 
TeV source HESS J1837— 069. At an assumed distance of 6.6 kpc by association with an adjacent 
massive star cluster, the efficiency of PSR J1838— 0655 converting spin-down luminosity to radiation is 
0.8% for the 2-10 keV ASCA flux, 9% for the 20-300 keV INTEGRAL flux and ~ 3% for the > 200 
GeV emission of HESS J1837— 069, making it a plausible power source for the latter. A Chandra X- 
ray observation resolves AX J1838. 0—0655 into a bright point source surrounded by a ss 2' diameter, 
centrally peaked nebula. The spectra of the pulsar and nebula are each well fitted by power laws, with 
photon indices r = 0.5(0.3 — 0.7) and T = 1.6(1.1 — 2.0), respectively. The 2—10 keV X-ray luminosities 
of the pulsar and nebula are Lpsr — 4.6 x 10 34 o% 6 ergs s _1 and Lpwn = 5.2 x 10 33 d 2 6 ergs s _1 . 
A second X-ray source adjacent to the TeV emission, AX J1837.3— 0652, is resolved into an apparent 
pulsar/PWN; it may also contribute to HESS J1837— 069. The star cluster RSGC1 may have given 
birth to one or both pulsars, while fueling TeV emission from the extended PWN with target photons 
for inverse Compton scattering. 

Subject headings: pulsars: individual (AX J1838. 0-0655, AX J1837.3-0652, PSR J1838-0655) - 
stars: neutron — supernova remnants — X-rays: stars 
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1. INTRODUCTION 

Surveys of the Galactic plane (Aharonian et al. 2006a) 
by the High Energy Stereoscopic System (H.E.S.S.) find 
that at least half of its sources (Funk 2007) can be iden- 
tified with supernova remnants (SNRs) or pulsar wind 
nebulae (PWNe). The sizes of those sources that are 
identified as PWNe are generally larger in TeV 7-rays 
than in X-rays, revealing a new view of the penetration of 
high-energy particles into the surrounding medium. The 
larger TeV nebulae are often displaced from the pulsar, 
possibly by the reverse shock of a supernova that ex- 
ploded in an inhomogeneous medium. They may also 
contain relic electrons from the more energetic, younger 
phase of the pulsar spin-down. The characteristic mis- 
match between X-ray and 7-ray sizes, together with the 
positional offsets of the X-ray and TeV emitting nebulae, 
are becoming familiar as new examples are found (Aha- 
ronian ct al. 2006b, 2006c, 2007). Understanding exactly 
how the spin-down luminosity of a particular pulsar pow- 
ers its TeV nebula requires time-dependent modeling and 
sometimes uncertain details of the local environment that 
must be determined at other wavelengths. The favored 
theoretical mechanisms, either inverse-Compton scatter- 
ing of ambient photons, or decay of neutral pions pro- 
duced in hadronic collisions of high-energy protons with 
a dense phase of the ISM, are being tested in several 
cases. For a recent review of these issues, see de Jager & 
Djannati-Atai (2008). 

A most instructive association is AX J1838. 0—0655 
with HESS J1837— 069. One of the first extended sources 
detected by H.E.S.S. (Aharonian et al. 2005, 2006a), it 
was considered unidentified until now even though X-ray 



studies of the field pointed to a coincident hard, steady 
X-ray source detected up to 300 keV by INTEGRAL 
(Malizia et al. 2005). AX J1838.0-0655 was seen by 
X-ray satellites spanning decades, beginning with Ein- 
stein (IE 1835.3-0658; Hertz & Grindlay 1988), and in- 
cluding ASCA (Bamba et al. 2003), BeppoSAX (Mal- 
izia et al. 2005), XMM-Newton (unpublished), and Swift 
(Landi et al. 2006), always with steady flux. Bamba 
et al. (2003) measured 1.1 x 10~ n ergs cm -2 s _1 from 
AX J1838.0-0655 in the 0.7-10 keV band, which is com- 
parable to the > 200 GeV flux from HESS J1837-069, 
3.3 x 10 -11 ergs cm -2 s _1 (Aharonian et al. 2006a). The 
absence of variability is an important clue to the nature 
of AX J1838. 0—0655, as is its hard spectrum. Malizia et 
al. (2005) fitted the combined INTEGRAL and ASCA 
spectrum with a power law of photon index T = 1.5, and 
conjectured that the source is a supernova product, more 
likely a PWN than a shell SNR. 

In §2 we describe our analysis of a Chandra im- 
age of HESS J1837-069 that resolves AX J1838.0-0655 
into a pulsar/PWN, and study its spectral proper- 
ties. We also address the nature of the other X-ray 
sources in the vicinity of HESS J1837— 069, one of which, 
AX J1837.3-0652, is evidently a second PWN. Moti- 
vated by the Chandra results, we proposed observations 
of AX J1838. 0-0655 with RXTE to search for its puta- 
tive pulsar. The resulting discovery of PSR J1838— 0655 
and measurement of its spin-down parameters are de- 
scribed in §3. We show that the rotational energy/spin- 
down luminosity of PSR J1838— 0655 are sufficient to 
power the TeV emission of HESS J1837-069, and we 
further argue in §4 that the adjacent massive star cluster 
RSGC1 (Figer et al. 2006; Davies et al. 2008) provides 
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Fig. 1. — Chandra ACIS-I broad-band X-ray image of HESS J1837— 069. Left: Image scaled to emphasize point sources listed in Table 1. 
Right: Image smoothed to emphasize diffuse emission. The ellipse represents the two-dimensional Gaussian Icr extent of the TeV emission 
(Aharonian et al. 2006a). The circle is the approximate extent of the red supergiant cluster RSGC1 (Figer ct al. 2006). The squares 
indicate the zoomed-in regions around AX J1838. 0—0655 and AX J1837.3— 0652 corresponding to the insets in Figure 2. 



ambient target photons for inverse Compton scattered 
TeV emission by PWN electrons from AX J1838.0-0655 
and possibly from AX J1837.3— 0652, whose pulsar may 
also have been born in the cluster. 

2. CHANDRA OBSERVATIONS AND RESULTS 

A 20 ks Chandra X-ray observation of 
HESS J1837-069 was made on 2006 August 19 
UT using the Advanced CCD Imaging Spectrometer 
(ACIS; Burke et al. 1997) operating in the full- frame 
TIMED/VFAINT exposure mode. This detector is 
sensitive to X-rays in the 0.3—12.0 keV energy range 
with a resolution of AE/E ~ 0.06 FWHM at 1 keV. 
The imaging system offers an on-axis spatial resolution 
of w 0.5", which is also the instrument pixel size, and 
increases to ~ 5" at 10' off-axis. A total of 19.9 ks 
of live-time was accumulated with a CCD frame time 
of 3.241 s (given the 1.3% readout deadtime). With a 
maximum count rate in a pixel of < 0.004 s _1 , photon 
pile-up can safely be ignored. As normal procedure, the 
spacecraft was dithered to average out the single-pixel 
spectral response. We used the standard processed and 
filtered event data. No time filtering was necessary 
as the background rate was stable over the course of 
the observation. All data reduction and analysis was 
performed using the CIAO (V3.4), FOOLS (V6.0.4), 
CALDB 3.4.2, and XSPEC (V12.2.1) X-ray analysis 
software packages. We followed the CIAO online science 
thread to create the exposure-corrected image. 

Figure 1 is the exposure-corrected map of the entire 
4-CCD ACIS-I field that is smoothed on two scales to 
highlight cither the point sources or the diffuse emission. 
There are 18 unambiguous point sources (listed in Ta- 
ble 1); two of these, sources 1 and 14, have associated 
X-ray nebulae. Since the image was centered on the TeV 
emission, AX J1838.0-0655 (source 14) fell 5.25' off-axis, 



near the gap between ACIS-I CCDs. Fortunately, expo- 
sure correction using the spacecraft dithering informa- 
tion fully recovers the fluxed image, allowing us to gen- 
erate an accurate spectrum and radial profile. 

2.1. AX J '1838. 0-0655 

Figure 1 shows that AX J1838.0-0655 falls along the 
long axis of the TeV extent 1 . Considering the off-axis lo- 
cation of the point source, we determine its precise loca- 
tion by fitting the off-axis point spread function following 
the CIAO recipe. Then, by registering other X-ray point 
sources with their optical or radio counterparts (see Ta- 
ble 1) we are able to refine the location of the point source 
in AX J1838.0-0655 to (J2000.0) R.A. = 18 h 38 m 03.13 s , 
deck = -06°55'33.4" with a la error radius of w 0.3". 
(A correction of —0.55" in R.A. and —0.09" in deck was 
applied.) 

Figure 2 shows an enlarged region of the Chandra im- 
age of AX J1838. 0—0655, and its radial profile, clearly 
resolving a nebula from a point source. For spectral 
analysis of AX J1838. 0—0655 we extract photons in the 
2—10 keV energy band, below which the spectra are 
highly attenuated by the large column density. For the 
point-source spectrum, we extract 2124 photons from a 
5" x 7" elliptical aperture centered on the source peak; 
the diffuse emission contributes a negligible background 
(21 counts) in this region. For the nebula we use a 1' 
radius extraction region, offset from the point source 
and centered at (J2000.0) R.A. = 18 h 38 m 03.1 s , deck = 
— 06°55'40" with a 7" x 9" region around the pulsar re- 
moved. The background for the nebula was extracted 

1 Figure 6 of Landi et al. (2006) that graphs the TeV extent on 
the Swift X-ray image has the wrong position angle for the ellipse, 
which doesn't do justice to the case for association. We reproduce 
here the orientation given by Aharonian et al. (2006a). 
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TABLE 1 

Point Sources in Chandra ACIS-I ObsID 6719 



No. 


X-ray 
R.A. (J2000) 


Position 
Decl. (J2000) 


Counts 


HR a 


Optical or 
R.A. (J2000) 


Radio Position 
Decl. (J2000) 


B2 b 
(mag) 


R2 b 
(mag) 


ID (Spectral Type) 


1 


lo 


QV Ol 91 
O / 


—06 52 30 6 


17 


+0.77 










AX 11 837 3—0652 

ilA -J lOO I . <J UUJi 


2 . 


lo 


Q7 99 1 A 
Ol ZZ.l^l 


— Ofi ^8 37 ^ 

UU UO O 1 .o 


83 


—0.36 


1 £ ^7 99 1 9fi 

lO Ol ^^.l^U 


— Ofi HQ 
uu uo oo.uy 


16.73 


14.30 




3 . 


1 8 
lo 


Q7 99 9fi 
O I ZZ.ZU 


— Ofi zL7 zLfl zL 

— UU rt 1 "rtU.rt 


59 


+0.82 












4 . . 


18 


37 37.37 


— 07 DO 09 fi 

— Uf UU UZr.U 


19 


—0.58 


1 R 17 ^7 ^74 

lO O 1 O ( .O ( 1 


— 07 DO 09 47 
U / UU uz.<± / 


16.50 


13.64 




5 . . 


18 


37 38.04 


— f)f> ^3 01 8 

— UU UO Ui.O 


25 


+0.92 












6 


18 


37 41.84 


— OtS 4fi 33 

UU '-tVJ J'J-U 


45 


+0.86 












1 .... 


18 


37 50.45 


-06 48 56.4 


85 


+0.98 


18 37 50.45 


-06 48 56.9 






GPSR5 25.320-0.098 


8.... 


18 


37 51.57 


-06 53 45.4 


274 


+0.99 


18 37 51.60 


-06 53 45.0 






GPSR5 25.252-0.139 


9.... 


18 


37 54.71 


-06 57 05.6 


23 


+0.91 












10. . . 


18 


37 58.00 


-06 53 31.6 


55 


+0.96 


18 37 57.99 


-06 53 31.0 






GPSR5 25.266-0.161 


11. .. 


18 


37 58.00 


-07 03 42.4 


44 


+0.85 












12. . . 


18 


37 58.72 


-06 48 22.5 


129 


-0.87 


18 37 58.760 


-06 48 22.39 


8.76 


7.75 


HD 171999 (G5) 


13. . . 


18 


37 59.97 


-06 49 22.4 


285 


+0.96 










Swift J1837.9-0649 


14. . . 


18 


38 03.13 


-06 55 33.4 


2114 


+0.92 










AX J1838.0-0655 


15... 


18 


38 10.59 


-06 51 31.9 


46 


+0.80 












16... 


18 


38 14.26 


-06 48 21.7 


240 


+0.65 












17... 


18 


38 15.27 


-06 47 58.9 


291 


+0.42 


18 38 15.231 


-06 47 58.78 




16.06 


W42 1 (05.5) 


18. . . 


18 


38 16.75 


-06 48 49.4 


126 


-0.06 











NOTE. — Includes all point sources that have signal-to-noise ratio > 3 in the 0.3—10 keV band, as determined by the CIAO 
software package source-detection tool wavdetect. Units of right ascension are hours, minutes, and seconds, and units of declination 
are degrees, arcminutes, and arcseconds. 



a Hardness ratio defined as HR = (N^ — N s ) / ( + N s ) , where N s and are the counts measured in the 0.3 — 2 keV and 2 — 10 kcV 
energy band, respectively. 
b Optical data from the USNO-B1.0 catalog (Monet et al. 2003). 




RADIAL DISTANCE (ARCSEC) RADIAL DISTANCE (ARCSEC) 

Fig. 2.— Chandra ACIS radial profiles of AX J1838.0-0655 (left) and AX J1837.3-0652 (right) compared to the scaled point spread 
function at their locations (FWHM ra 2"; solid line). Surrounding each point source, a clear flux excess is evident over the field background 
flux level (horizontal dashed line) torss 60" in the case of AX J1838. 0—0655. For AX J1837.3— 0652, a larger but lower surface brightness 
extension is seen. Insets: Zoom-in on the Chandra images. The intensity is scaled to highlight the diffuse emission; the point sources are 
saturated. 



from a concentric annulus of radii 1' < r < 2'. The neb- 
ula so extracted comprises 483 photons after subtracting 
background (33%) and excising the point source. 

Spectra from the point source and nebula region were 
grouped with a minimum of 50 counts per spectral chan- 
nel and fitted using XSPEC. The two spectra were fitted 
simultaneously using absorbed non-thermal power laws 
with their hydrogen column density parameters linked; 
the results are summarized in Figure 3. The chosen 
model produced excellent fits (x 2 = 0.9 for 57 DoF) 
for a common column density of ATjj ~ 4.5(3.7 — 5.2) x 
10 22 cm~ 2 (90% confidence interval is used throughout). 



The spectrum of the putative pulsar has a photon index 
of T = 0.5(0.3 — 0.7), similar to that found for the pulsar 
PSR J1811-1926 in the young SNR Gil. 2-0. 3, and at 
the low (hard) end of that expected for a high-i? rotation- 
powered pulsar (> 4 x 10 36 erg s" 1 ; cf. Gotthelf 2003). 
For the PWN emission we measure a photon index of 
T = 1.6(1.1-2.0), which is also at the low end for a PWN. 
The unabsorbed 2—10 kcV fluxes for the putative PSR 
and PWN are F PSR = 8.8 x 10~ 12 ergs cur 2 s^ 1 and 
Fpwn = 1-0 x 10~ 12 ergs cm -2 s _1 , respectively. These 
results are consistent with the ASCA composite spec- 
trum, having photon index T = 0.8, A^h = 4 x 10 22 cm+ 2 , 
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Fig. 3.— Chandra ACIS spectra of PSR J1838-0655 and its 
nebula, comprising AX J1838. 0—0655, both fitted to an absorbed 
power-law model (see text). Residuals for each spectrum are shown 
in units of standard deviation. 



and F(0.7- 10 kcV) = 1.1 x 1CT 11 ergs cmT 2 s" 1 (Bamba 
et al. 2003). The flux ratio F PWN /F PSR = 0.11 is at the 
low end of the distribution of PWNe observed by Chan- 
dra, which range from 0.1—10 as compiled by Kargaltsev 
& Pavlov (2008). 

The Chandra position of the point source rules out the 
optical/IR candidate suggested by Landi et al. (2006) on 
the basis of a less accurate Swift position. We deter- 
mine an optical upper limit of R > 23.0 from an image 
obtained on the 2.4m Hiltner telescope of the MDM Ob- 
servatory (Fig. 4). One does not expect to detect a pulsar 
having the large distance and extinction implied by the 
fitted Nh in an optical image of this depth. 

2.2. AX J 1837. 3- 0652 

In addition to AX J1838. 0—0655, a fainter point source 
(number 1 in Table 1) of 17 photons is surrounded 
by a large, diffuse nebula. It is evidently also a 
pulsar/PWN system. This source was detected by 
ASCA (AX J1837.3-0652; Bamba et al. 2003), with 
a 0.7-10 kcV flux of 1.8 x 10~ 12 ergs cm~ 2 s"\ and 
by Swift (sources 6 and 7 of Landi et al. 2006). 
It is difficult to measure here due to its low surface 
brightness. It is more extended than AX J1838. 0—0655 
(see Figs. 1 and 2) but has comparable total nebu- 
lar flux, w 1.2 x 10~ 12 ergs cm~ 2 s _1 (unabsorbed) 
in the 2-10 keV band. Unlike AX J1838.0-0655, in 
AX J1837.3-0652 the PWN flux dominates over the 
(putative) pulsar point source. Its spectral parameters 
are largely unconstrained, with 0.7 < T < 3.6 and 
N H = (2 - 12) x 10 22 cm" 2 . If AX J1838.0-0655 is 
to be advanced as the origin of HESS J1837— 069, then 
it cannot be ruled out that the PWN AX J1837.3-0652 
makes a significant contribution. 

2.3. Additional Chandra X-ray Sources 



Sources 16—18 of Table 1 lie in the northeast corner 
of the Chandra image, and are evidently associated with 
the W42 star forming region, with source 17 being the 
central 05.5 star in the cluster (Blum et al. 2000). Three 
additional sources (numbers 2, 4, and 12) are also stars, 
which is consistent with their softer X-ray spectra. HD 
171999 in particular is a high-proper-motion star. Its op- 
tical position listed in Table 1 is calculated for the epoch 
of the Chandra observation. None of the other sources 
have optical/IR counterparts in the Digitized Sky Survey 
or 2MASS; they are probably active galactic nuclei. 

Positions for radio sources denoted by GPSR5 desig- 
nations in Table 1 are measured from the 1.4 GHz VLA 
Multi-Array Galactic Plane Survey (White et al. 2005; 
Helfand et al. 2006) 2 . The two radio sources nearest 
the TeV emission (numbers 8 and 10) were discussed by 
Figer et al. (2006) and Trejo & Rodriguez (2006). The 
latter authors proved that they are extragalactic based 
on the kinematics of their H I absorption spectra. Their 
X-ray positions coincide with compact components of 
the radio sources, presumably their host galaxy nuclei 
or quasars. We have no reason to hypothesize that any 
of these sources are related to HESS J1837— 069. 

3. RXTE OBSERVATIONS AND RESULTS 

The point-source flux measured by Chandra in 
AX J1838. 0—0655 is sufficient that a pulsar search with 
RXTE (Bradt et al. 1993) was deemed feasible. Ac- 
cordingly, we requested six observations in total, pointed 
at AX J1838.0-0655, that spanned 2008 February 17 - 
March 5 UT to search for the expected pulsar signal and 
develop a phase-connected timing solution. A log of the 
observations is given in Table 2. The data used here were 




Fig. 4. — A 600 s _R-band exposure of the center of 
AX J1838. 0-0655 obtained on the 2.4m Hiltner telescope on 2005 
July 3. The displayed field is 30" X 30". North is up, and east is 
to the left. The Chandra position of PSR J1838— 0655 is indicated 
by a circle of radius 0.7", approximately its 2 <r uncertainty. The 
magnitude limit is R = 23.0. 

2 See http://third.ucllnl.org/gps/. 



Discovery of a Pulsar Associated with HESS J1837-069 



5 



TABLE 2 

Log of RXTE Observations and Period Measurements 



TABLE 3 
Timing Parameters of PSR J 1838- 



0655 



Date 


Exp. /Span 


Start Epoch 


Period a 


zl 


CUT) 


(ks) 


(MJD) 


(ms) 


2008 Feb 17 


9.1/14.0 


54513.825 


70.49820(3) 


277 


2008 Feb 19 


20.1/35.0 


54515.476 


70.498218(6) 


726 


2008 Feb 22 


5.4/8.0 


54518.609 


70.49828(8) 


153 


2008 Feb 26 


5.5/8.2 


54522.004 


70.49823(6) 


186 


2008 Mar 2 


5.3/8.0 


54527.378 


70.49829(6) 


223 


2008 Mar 5 


6.3/8.6 


54530.185 


70.49830(6) 


206 



a Period derived from a Z\ test. The Leahy et al. (1983) uncertainty 
on the last digit is in parentheses. 



collected with the Proportional Counter Array (PCA; 
Jahoda et al. 1996) in the GoodXenon mode with an 
average of 3.6 out of the five proportional counter units 
(PCUs) active. In this mode, photons are time-tagged to 
0.9 /zs and have an absolute uncertainty better than 100 
[is (Rots et al. 1998). The effective area of five combined 
detectors is about 6500 cm 2 at 10 keV with a roughly 
circular field-of-view of ~ 1° FWHM. Spectral informa- 
tion is available in the 2—60 keV energy band with a 
resolution of ~ 16% at 6 keV. 

Standard time filters were applied to the realtime PCA 
data, which rejects intervals of South Atlantic Anomaly 
passages, Earth occultations, and other periods of high 
particle activity. Our data contained several short (<~ 
200 s) intervals of detector arcing which were excised. 
The realtime data required manually adjusting the good 
data intervals to properly correspond to actual data col- 
lection times for each PCU. Table 2 lists the net expo- 
sure times and total spans of each data set. The pho- 
ton arrival times were transformed to the solar system 
barycenter in Barycentric Dynamical Time (TDB) using 
the JPL DE200 cphemcris and the Chandra measured 
coordinates. 

3.1. Timing Analysis 

We restricted the analyzed data to the 2—20 keV en- 
ergy range (PCA channels 2—50) from the top Xenon 
layer of each PCU to optimize the signal-to-noise. A 
Fast Fourier Transform (FFT) of the first observation 
(February 17) revealed a significant signal of power 
S = 117 at a period of P = 70.498 ms 3 . For 
the 5.7 x 10 6 search elements this corresponds to es- 
sentially nil false detection probability. There is no 
doubt that PSR J1838— 0655 is the compact object 
in AX J1838. 0—0655, given the morphological evidence 
from Chandra and because AX J1838. 0—0655 is the 
brightest known X-ray source in the RXTE field of view. 

Based on the strength of the signal in the February 17 
and 19 detections of PSR J1838— 0655, the next four ob- 
servations were scheduled from February 22 - March 5 to 
develop a fully phase-coherent timing solution including 
the previous data sets. For each epoch, we extracted the 
pulse profile corresponding to the peak period as deter- 

3 While planning these observations, the signal was first found 
by C. Markwardt (personal communication) in summed FFTs of 
hundreds of short, serendipitous RXTE dwells on the source taken 
during 2004-2008. However, the individual ~ 100 s exposures did 
not contain sufficient statistics to accurately measure the period or 
derive an ephemeris. 



Parameter Value 

R.A. (J2000) a 18 h 38 m 03.13 s 

Decl. (J2000) a -06°55'33.4" 

Epoch (MJD TDB) 54522.00000012 

Period, P (ms) 70.498243969(54) 

Period derivative, P 4.925(29) X 10~ 14 

Range of timing solution (MJD) 54513-54530 

Characteristic age, t c (kyr) 22.7 

Spin-down luminosity, E (ergss -1 ) . . 5.5 X 10 36 

Surface dipolc magnetic field, B s (G) 1.9 X 10 12 

NOTE. — TEMPO la uncertainties given in parentheses. 
a Chandra ACIS-I position from Table 1. 



5 10 15 

TIME (DAYS SINCE MJD 54513) 

Fig. 5. — Phase residuals of the six observations 
PSR J1838— 0655 from the quadratic cphemcris of Table 3. 



of 



mined by the Z\ test (Buccheri et al. 1983). The result- 
ing profiles were cross correlated, shifted, and summed 
to generate a master pulse profile template. Individual 
profiles were then cross correlated with the template to 
determine the time of arrival (TOA) and its uncertainty 
in phase at each epoch. These TO As were iteratively fit- 
ted to a quadratic ephemeris using the TEMPO software. 
We started by fitting just TOAs from the first three 
epochs, which straddle the longest observation (Febru- 
ary 19), to a linear solution, and then added the last 
three TOAs to the fit one at a time. At each step we 
found that the new TOA would match to < 0.1 cycles 
the predicted phase derived from the previous set. The 
complete resulting ephemeris is presented in Table 3, and 
the phase residuals are shown in Figure 5. The residuals 
are all less than 0.02 cycles and appear to be random 
within their statistical uncertainties. Thus, inclusion of 
higher order terms is not indicated. An independent test 
of the uniqueness of the timing solution was performed 
using a Z\ search on a two-dimensional grid of P and P, 
from which we obtained consistent ephemeris parameters 
and errors. 

Figure 6 displays the summed pulse profile using all 
the 2—20 keV data folded on the final cphemcris. It 
has a symmetric, blended double-peaked structure. The 
measured pulsed fraction of 4.8% indicates a intrinsic 
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Fig. 6.— RXTE folded light curve of PSR J1838-0655 in 
the 2—20 keV band. Phase zero corresponds to the epoch of the 
cphcmcris in Table 3. 



modulation of ^ 50% after allowing for the estimated 
internal and Galactic background in the PCA. We see no 
energy dependence of the pulse profile when subdividing 
the 2-20 keV band. 

The spin-down parameters of PSR J1838— 0655, partic- 
ularly its E — 5.5 x 10 36 ergs s _1 , are in accord with ex- 
pectation from the Chandra results on AX J1838. 0—0655 
described in §2. Assuming a distance of 6.6 kpc as dis- 
cussed in §4 the ratio L x (2-10 keV)/£ = 8xl0~ 3 for the 
pulsar, and 1 x 10~ 3 for its PWN. The 20-300 keV flux of 
AX J1838.0-0655 measured by INTEGRAL (9 x 10" 11 
ergs cm~ 2 s _1 ; Malizia et al. 2005) corresponds to a lu- 
minosity of 0.09 Edl 6 , and the > 200 GeV flux from 
HESS J1837-069 (3.3 x 10" 11 ergs cm" 2 s" 1 ; Aharonian 
et al. 2006a) corresponds to 0.03 Edl 6 . 

3.2. RXTE Spectral Analysis 

The spectrum of the pulsed flux from PSR J1838-0655 
can be isolated using phase-resolved spectroscopy. We 
used the fasebin software to construct phase-dependent 
spectra based on the ephemeris of Table 3. For each 
epoch and PCU we constructed spectra from the top 
Xenon layer only and combined them to produce sin- 
gle spectra per PCU for the entire set of observations. 
Similarly, standard PCA responses for each PCU were 
generated at each epoch and averaged. In fitting the 
pulsed flux, the unpulsed emission provides a near per- 
fect background estimate. The spectra were divided into 
two groups of 0.5 cycles each to represent the "off-peak" 
and "on-peak" emission and were fitted in the 2 — 20 keV 
range using XSPEC. The spectra from the four active 
PCUs were analyzed, both separately and combined into 
a single spectrum, which gave similar results. 

A simple absorbed power-law model was used with the 
interstellar absorption held fixed at Ah = 4.5 x 10 22 cm~ 2 
determined from the Chandra fit; leaving the Nn uncon- 
strained results in a larger uncertainty in T. The result- 
ing best-fit photon index is T = 1.2(1.1 — 1.3). The unab- 
sorbed 2—10 keV pulsed flux is 9 x 10~ 12 ergs cm~ 2 s _1 , 
which represents, within statistics, all of the point- 
source flux of PSR J1838— 0655 measured by Chandra 
(8.8 x 10~ 12 ergs cm~ 2 s _1 ), indicating that its intrin- 
sic pulsed fraction is ~ 100%. The corresponding pulsed 
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Fig. 7.— RXTE spectrum of pulsed flux from PSR J1838-0655 
obtained by subtracting the off-peak spectrum from the on-peak 
spectrum, and fitting to an absorbed power-law model (see text). 



luminosity (assumed isotropic) is 4.6 x 10 34 d\ 6 ergs s" 1 . 

The RXTE power-law slope derived from the pulsed 
emission is slighlty flatter than T = 1.5 reported by Mal- 
izia et al. (2005) for the joint fit to ASCA and INTE- 
GRAL spectra of AX J1838.0-0655. We also find that 
this result is at odds with the Chandra spectra presented 
in §2.1, which indicate T — 0.5 ± 0.2 for the dominant 
pulsar component. The difference is sufficiently signifi- 
cant to suggest that a steepening of the spectrum of the 
pulsar occurs in the 8—15 keV band. 

4. DISCUSSION 

Graphic evidence of the relationship of 
AX J1838.0-0655 to HESS J1837-069 is provided 
by the Chandra ACIS-I image, which resolves 
AX J1838. 0—0655 into a point source plus a cen- 
trally peaked 2' diameter nebula near one end of the 
major axis of the TeV source. Interestingly, a second 
pulsar/PWN system is evident in the Chandra image 
corresponding to a weaker source previously detected 
as AX J1837.3— 0652. Although AX J1837. 3-0652 
contains a point-source component ^100 times weaker 
than that of AX J1838. 0-0655, the nebular fluxes of 
the two sources are virtually identical. This is not 
too surprising, as the ratio of pulsar X-ray luminosity 
to spin-down power is such a widely ranging quantity 
(Kargaltsev & Pavlov 2008). The mere presence of the 
second PWN is of relevance here because it indicates 
another likely energetic pulsar of E > 4 x 10 36 ergs s _1 
(Gotthelf 2003), which implies that the two sources are a 
priori comparably suited to powering the TeV emission. 
Since AX J1837.3— 0652 is located near the opposite end 
of the elongated TeV source from AX J1838.0-0655, the 
possibility that HESS J1837-069 comprises two PWNe, 
perhaps born in the same star cluster, is an intriguing 
question for further investigation. 

The location of AX J1838. 0—0655 near an unusual 
cluster of red supcrgiant stars (RSGC1; Figer et al. 2006) 
suggests a distance of 6.6 kpc that Davies et al. (2008) 
determined from the radial velocity of the cluster. Davies 
et al. also measured the infrared extinction to the clus- 
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ter, Ak = 2.6 mag, which is equivalent to Ay — 23.5, 
or TVh = 4 x 10 22 cm' 2 according to the conversion 
Nn/A v = 1.79 x 10 21 cm- 2 mag" 1 (Predehl & Schmitt 
1995). Our X-ray fitted Ah is consistent with this. It is 
plausible that PSR J1838— 0655 was born in this young 
cluster, which is centered 3' north-northwest of the X- 
ray source (see Figure 1). At this distance, a 23 kyr old 
neutron star born in the center of the cluster requires 
a transverse velocity of 245 km s _1 to reach its present 
location, less if it was born in the outskirts. 

It is possible that the TeV PWN is displaced in one 
direction from the pulsar by an asymmetric reverse 
shock resulting from a supernova that exploded ini- 
tially in an inhomogeneous medium, an explanation sug- 
gested by Blondin et al. (2001) for the similar appear- 
ance of the Vela X remnant (Aharonian et al. 2006b). 
Other likely associations that show such displacements 
are PSR B1823-13/HESS J1825-137 (Aharonian et 
al. 2006c; Pavlov et al. 2008), PSR J1809-1917/HESS 
J1809-193 (Aharonian et al. 2007; Kargaltsev & Pavlov 
2007), PSR J1718-3825/HESS J1718-385 (Aharonian 
et al. 2007; Hinton et al. 2007), and possibly PSR 
J1617-5055/HESS J1616-508 (Aharonian et al. 2006a; 
Landi et al. 2007). 

Based on the number of red supergiants in the clus- 
ter RSGC1, Davies et al. (2008) estimate that it has 
a total mass of (2 — 4) x 10 4 M Q , the highest of any 
young cluster in the Galaxy. Figer et al. (2006) calcu- 
lated that it should produce a supernova every 40—80 kyr 
on average, which is frequent enough statistically to ex- 
pect a pulsar of the 23 kyr characteristic age that we 
measured for AX J1838. 0—0655. In fact, it is statis- 
tically suggestive that the previous pulsar born in this 
cluster is AX J1837.3— 0652, and that it contributes to 
HESS J1837-069. Theoretically, TeV sources powered 
by pulsars can live for up to 100 kyr as inverse Comp- 
ton emitters in regions of low (/xG) -B-fields (de Jager & 
Djannati-Atai 2008). 

A plausible mechanism of the TeV emission from 
HESS J1837— 069 is inverse Compton scattering of am- 
bient photons off relativistic electrons in an extended 
PWN, with the target IR/optical photons supplied by 
the star cluster. Considering the luminosity of only the 
IR detected supergiants of RSGC1 as a lower limit on 
the seed photons (« 2.8 x 10 6 L Q : Figer et al. 2006; 
Davies et al. 2008), at a distance of 3 x 10 19 cm from 
the cluster the energy density of optical/IR photons in 
HESS J1837-069 is - 20 eV cm" 3 . If the pulsar injected 
10 37 e s _1 into this region for the past 10 4 yr with av- 
erage 7 = 10 6 , then the inverse Compton luminosity can 
be < 2.7 x 10 36 ergs s _1 . This approximation does not 
account for the Klein-Nishina suppression of the cross 
section, which would reduce it. Nevertheless, it appears 
sufficient to supply the 1.7 x 10 35 dg 6 ergs s" 1 luminos- 
ity of HESS J 1837— 069, as it does not even invoke the 
higher luminosity from the pulsar in its younger phase, 
which energy may still be present as relic electrons in the 
extended TeV nebula. 



5. CONCLUSIONS 

Using Chandra, we resolved the X-ray image and 
spectrum of AX J1838. 0—0655 into its pulsar and 
PWN components, with 2—10 kcV luminosities of 
Lpsr = 4.6 x 10 34 dg 6 ergs s _1 and L PW n = 5.2 x 
10 33 d\ 6 ergs s _1 . Led by this evidence we searched for 
pulsations using RXTE and discovered PSR J1838-0655 
in AX J1838.0-0655 at a period of 70.5 ms. Timing 
its spin-down, we found E — 5.5 x 10 36 ergs s" 1 , a 
value that is consistent with pulsars of similar X-ray 
luminosity and X-ray nebular properties. The 20—300 
kcV luminosity observed by INTEGRAL (Malizia et al. 
2005) is w 4.7 x 10 35 d\ 6 ergs s" 1 , and is likely to be an 
extension of the hard pulsar spectrum rather than the 
PWN. The 0.2-20 TeV luminosity of HESS J1837-069 
is 1.7 x 10 35 <i§ 6 ergs s _1 , which is 3% of the pulsar E, 
while the TeV emission is several times larger in extent 
than the X-ray PWN. PSR J1838-0655 is located at one 
end of the elongated TeV source. 

The 23 kyr characteristic age of PSR J1838-0655 is 
comparable to those of the middle-aged pulsars Vela and 
PSR B1823— 13, with which it shares the characteristic 
of an X-ray PWN displaced from the center of a larger 
TeV nebula. Such pulsars are old enough that the dis- 
placement can be enforced by a reverse supernova shock 
propagating in an inhomogeneous interstellar medium. 
Alternatively or in addition, we may be seeing relic elec- 
trons from asymmetric diffusion out of the PWN in an 
earlier phase. 

We also discovered a Chandra point source inside a sec- 
ond diffuse X-ray source adjacent to HESS J1837— 069. 
Probably an older pulsar/PWN, AX J1837.3-0652 may 
also contribute to HESS J1837— 069, explaining its elon- 
gated shape as a blended source. This scenario may 
be testable with spatially resolved spectroscopy of the 
TeV emission, especially if the age of the putative 
second pulsar is significantly different from the first. 
Although not yet conclusively identified as a pulsar, 
AX J1837.3— 0652 should have similar spin-down lumi- 
nosity as AX J1838. 0—0655. A deep search for radio 
pulsations from AX J1837.3— 0652 is warranted to de- 
termine its spin-down parameters, it being too faint for 
an X-ray search. 
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